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We search for resonant WW or W Z production using up to 5.4 fb" 1 of integrated luminosity 
collected by the DO experiment in Run II of the Fermilab Tevatron Collider. The data are consistent 
with the standard model background expectation, and we set limits on a resonance mass using 
the sequential standard model (SSM) W' boson and the Randall-Sundrum model graviton G as 
benchmarks. We exclude an SSM W' boson in the mass range 180 — 690 GeV and a Randall- 
Sundrum graviton in the range 300 — 754 GeV at 95% CL. 

PACS numbers: 12.60.Cn, 13.85.Rm, 14.70.Kw, 14.70.Pw 



The standard model of particle physics is expected to 
be a low energy effective theory valid for particle interac- 
tions below the TeV scale. Above this scale, extensions 
to the standard model (SM) augment the existing parti- 
cle content, leading to enhanced production of many final 
states at colliders. Specifically, the production and de- 
cay of massive charged or neutral particles can produce 
an excess of W boson pairs for neutral particles or W and 
Z boson pairs for charged particles. 

In this Letter, we search for resonant 
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WW and WZ production using data collected by 
the DO detector from 1.96 TeV pp collisions produced 
by the Fermilab Tevatron Collider. We use a sequential 
standard model (SSM) 0, Q W' boson as bench- 
mark for a WZ resonance and a Randall-Sundrum 
(RS) 3 5] graviton (G) resonance for the WW final 
state. 

There are two recent direct searches for 
WZ or WW resonances by the CDF and DO col- 
laborations that exclude WZ resonances with mass 
below 516 and 520 GeV, respectively, and an RS 
graviton G — > WW resonance with mass less than 
607 GeV Indirect searches for new physics in the 

WW and WZ diboson systems through measurements 
of the triple gauge couplings also show no deviation 
from the SM predictions [8h10|. The DO collabora- 
tion also excludes M(W') < 1.00 TeV [0], when 
assuming the W' boson decays as in the SM, and 
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M{G) < 1.05 TeV [l2j when assuming G decays into 
77 or ee. The search for these resonances in the diboson 
decay channel covers the possibility that their coupling 
to leptons may be lower than the value predicted by the 
SM. 

We obtain a combined result based on three in- 
dependent searches: two new searches for resonant 
WW/WZ production with at least one jet and exactly 
one or two leptons in the final state using 5.4 fb" 1 of 
integrated luminosity and one search previously done on 
4.1 fb -1 of integrated luminosity with three leptons in the 
final state @. We use data collected by the DO experi- 
ment. A detailed description of the DO detector can be 
found in [l3j ; we only give a brief description here. The 
innermost region is the tracking detector, which consists 
of silicon microstrip and central fiber trackers, both of 
which are surrounded by a solenoidal magnet produc- 
ing a 2 T magnetic field. Charged particle tracks are 
formed from signals in these detectors. Surrounding the 
tracking detector are electromagnetic (EM) and hadronic 
calorimeters, both of which use liquid argon as the active 
medium. The calorimeters are housed in three cryostats 
that define the central region as < 1.1 [l4[ and two 
endcap regions as 1.5 < \r)\ < 4. Electrons are recon- 
structed in the EM calorimeter as isolated energy clus- 
ters, matched to tracks, and with a shower shape that is 
consistent with that of an electron. Jets arc also formed 
in the calorimeters as clusters of energy in a cone with 
radius 1Z = 0.5 [15j |. Finally, surrounding the calorime- 
ters are central and forward muon systems in three layers 
consisting of precision wire chambers and fast scintilla- 
tors used for triggering. Coverage of the muon system 
extends to |?/| « 2. Located between the first and second 
layer of the muon system is a 1.8 T toroidal magnet which 
allows an independent muon momentum measurement . 
A muon candidate is reconstructed as the combination 
of tracks in the muon system and the inner tracking de- 
tector and is required to be isolated from other tracks or 
calorimeter energy deposits. 

We employ a Monte Carlo (MC) simulation to model 
all background processes except backgrounds from events 
not involving the decay of a W or Z boson, such as 
multijet production. The SSM W' boson, RS graviton, 
and the SM diboson processes are simulated using the 
PYTHIA [l6| event generator, which generates the tree- 
level matrix element process and simulates subsequent 
particle showering and hadronization effects. Back- 
grounds from ti, W+jets with W — > iv, and Z+jets 
with Z — !> ££ are modeled using the ALPGEN [l7| gen- 
erator and single top quark production is modeled with 
the COMPHEP [l8| generator. All generators are inter- 
faced with pythia for showering and hadronization. The 
ALPGEN-generated samples make use of the MLM [l!| jet- 
parton matching scheme to improve the jet multiplicity 
modeling. All MC samples are passed through a GEANT- 
based [20| simulation of the DO detector and overlaycd 



with data events from randomly selected bunch cross- 
ings to simulate multiple pp interactions within a single 
event. The signal samples are generated in exclusive fi- 
nal states with diboson resonance masses between 180 
and 1250 GeV in 10 GeV steps up to 200 GeV and then 
50 GeV steps above 200 GeV, using the CTEQ6L1 [HI 
parton distribution functions. No interference between 
the SM W boson and the SSM W boson production 
is included in the simulation since the effect is negligi- 
ble [22|] • All MC samples are normalized such that the 
predicted yield is equal to the production cross section 
multiplied by the integrated luminosity of data. The 
W+jets and Z+jets samples are scaled to the product 
of the cross section calculated by ALPGEN and the k- 
factor defined as the ratio of the next-to-leading order 
(NLO) and leading order (LO) cross sections, which is 
computed by MCFM [231 ] . The ti events are normal- 
ized to a next-to-next-to-leading order (NNLO) calcula- 
tion [24[ with m t = 172.5 GeV. Finally, the diboson sam- 
ples are normalized to the NLO cross section predicted 
by MCFM, and the signal W' boson samples are normal- 
ized to the NNLO cross sections [25j. The RS graviton 



samples are normalized to the pythia- level cross section 
multiplied by a fc-factor of 1.3 [2(|. 

Events in this search are placed in three mutually 
exclusive categories, thus maximizing signal sensitivity 
to each WW and WZ decay channel. The first cate- 
gory contains events with a leptonic decay of the W bo- 
son and hadronic decay of the W or Z boson. Events 
must contain exactly one electron or muon with trans- 
verse momentum pt > 20 GeV, cither one or two jets 
with px > 20 GeV, and missing transverse energy Jj] T > 
20 GeV [27J]. These events were collected using triggers 
that require the presence of a high px lepton. Events with 
charged leptonic decays of the Z boson and hadronic de- 
cays of the W boson comprise the second category. These 
events must contain exactly two electrons or muons and 
exactly one or two jets with the same px thresholds as 
the first category. We require $ T < 50 GeV to remove 
mismeasured events and a dilepton mass between 70 and 
110 GeV to select Z boson events. Both single lepton 
and dilepton triggers were used to collect events in this 
category. Fully leptonic decays of the WZ system con- 
stitute the final selection category. In this category any 
combination of three leptons (eee, /i/i/i, e/z/z, ee//) with 
Pt > 20 GeV for each lepton is accepted. Additionally, 
]JJt > 30 GeV is required. Events in this channel were 
collected using the same set of triggers as the dilepton 
channel. More details of the trilepton analysis have been 
presented in a previous Letter [6[ . 

In the first two selection categories the background 
after the initial event selection is dominated by W or 
Z boson+jets, followed by multijet, ti, single top quark, 
and diboson production. The multijet background, in 
both single lepton and dilepton events, is modeled us- 
ing data that fail the final lepton quality selection crite- 
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ria. In single lepton events, the relative fraction of mul- 
tijet (/afce-lcpton) background and all other backgrounds 
(real- lepton) is determined by measuring the relative 
rates at which each background type satisfies two differ- 
ent lepton quality criteria. A sample of Z — > it events is 
used to measure the reaZ-lepton efficiency, and a sample 
of events without significant missing transverse energy 
($ T < 20 GeV) is used to measure the /afce-lepton ef- 
ficiency. In dilepton events, the sum of the MC-based 
backgrounds and the multijet background is normalized 
to data in the dilepton mass region between 40 and 70 
GeV. The signal acceptance, in both single lepton and 
dilepton events, is estimated from the MC and corrected 
for all data/MC differences and the estimated trigger se- 
lection efficiency. 

The W and Z bosons result from the decay of a mas- 
sive resonance and are therefore highly boosted. We ex- 
ploit this property in the single lepton channel by re- 
quiring that the pr of the lepton- $ T system be greater 
than 100 GeV and the azimuthal angle between the lep- 
ton and ifjrp be less than 1.5 radians. In the dilepton 
channel, we require that the dilepton pair pt be greater 
than 100 GeV. The angular distance A1Z between the 
two leptons must be less than 1.5, where the angular dis- 
tance A1Z = y/ (Aif) 2 ± (A(j>) 2 is defined as the distance 
between the two objects in the pseudorapidity and az- 
imuthal angle plane. The hadronically decaying W or 
Z boson in both single lepton and dilepton events will 
also be highly boosted. Given the extended size of jets, 
the two jets from the hadronic decay of W or Z bosons 
with sufficient transverse momentum may be merged in a 
single jet whose mass, defined in terms of the jet energy 
Ej and momentum pj as rrij = \J E 2 — |pj| 2 , corresponds 

to the original boson mass. Thus, we further select events 
having a single jet with rrij greater than 60(70) GeV when 
searching for a hadronically decaying W(Z) boson. If no 
jets satisfy these requirements, two isolated jets must be 
reconstructed in the event with a dijet mass between 60 
and 105 GeV (70 and 115 GeV) for a W{Z) boson decay. 
In each channel, the AH between the jets must be less 
than 1.5 radians. 

We increase the search sensitivity by subdividing all 
search channels into "high" -mass and "low" -mass sig- 
nal regions, where the mass refers to the assumed sig- 
nal mass (M rcs ) and high-mass and low- mass are defined 
as M res > 450 GeV and < 450 GeV, respectively. The 
low-mass signal region is composed of all events that 
satisfy the single lepton and dilepton selection require- 
ments. In the high-mass signal region we additionally 
require the lepton- $ T system px > 150 GeV and the az- 
imuthal angle difference between the lepton and $ T to be 
less than 1.0 radian for single lepton events. The dilep- 
ton high-mass selection requires the dilepton px to be 
greater than 150 GeV and the ATI between the two lep- 
tons must to be less than 1.0 radian. Tabic Q] displays 



TABLE I: Expected event yields for the high-mass single lep- 
ton and dilepton selection samples. Uncertainties are com- 
bined statistical and systematic uncertainties on the back- 
ground yields. 



Process 


Single lepton sample 


Dilepton sample 


Z+jets 


3.6 ±0.2 


7.9 ±0.8 


W+jets 


124.5 ±20.3 


< 0.01 


Top 


22.9 ±2.5 


< 0.01 


Multijet 


4.6 ±0.3 


< 0.01 


Diboson 


27.6 ± 1.4 


0.8±0.1 


Background sum 


183.2 ±24.5 


8.7 ±0.8 


Data 


174 
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the estimated background yields, the expected numbers 
of signal events, and the numbers of observed data events 
after the high-mass selection. Figs. Q] and [2] compare the 
data with the estimated backgrounds in the single lepton 
and dilepton channels using the reconstructed resonance 
mass and transverse mass [28[ , respectively. In both Fig- 
ures, the absence of events below 350 GeV is a result of 
the high-mass event selection. 

The dominant systematic uncertainties on the back- 
ground normalization and signal acceptance in the dilep- 
ton channel are mostly due to Z±jets modeling (10%) 
and to jet energy resolution (3%) effects. The main un- 
certainties in single lepton events are W±jets modeling 
(15%), the tt cross section (10%), and the integrated lu- 
minosity (6.1%) [H]. In single jet events, the principal 
uncertainty is the jet-mass modeling. To determine it, a 
control region is defined from events satisfying the initial 
event selection, but failing the lepton- or the dilepton 
system px > 100 GeV requirement. The relative differ- 
ence between the background prediction and the data is 
10% for jet masses below 30 GeV and rises to 25% for 
masses near 60 GeV. We do not observe any event with 
jets having a mass greater than 60 GeV in the control 
region. 

No statistically significant excess of the data over the 
background prediction is observed. Thus, we set limits 
on the production cross section multiplied by the branch- 
ing ratio using a modified frequentist approach [3(|- In 
this method a log-likelihood ratio (LLR) test statistic [3l| 
is formed using the Poisson probabilities for estimated 
background yields, the signal acceptance, and the ob- 
served number of events for all resonant mass hypothe- 
ses. Confidence levels are derived by integrating the 
LLR in pseudo-experiments using both the signal plus 
background hypotheses {CL s+ i,) as well as the back- 
ground only hypothesis (CL&). In the modified frequen- 
tist approach, the excluded production cross section is 
computed as the cross section for which CL S , defined as 
CL s+ b/CLb, is equal to 0.05. Limits on the WZ reso- 
nance are set using the distribution of the reconstructed 
mass in dilepton events and of the reconstructed trans- 
verse mass in single lepton events, using a bin size of 50 
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D0, 5.4 fb _1 




• Data 

□ W+jets 

□tt 

□ ww+wz 

— SSM W'(600) 
■ ■ RS G(600) 



D0, 4.1-5.4 fb" 1 



300 400 500 600 700 800 
WW or WZ Transverse Mass [GeV] 

FIG. 1: (color online) The reconstructed WW or WZ trans- 
verse mass distribution for the data (points) with statistical 
error bars and the estimated W+jets, tt, and diboson back- 
grounds in the single lepton channel. The predicted SSM 
W' distribution for M(W') = 600 GeV and the RS gravi- 
ton distribution for M(G) = 600 GeV are shown as solid and 
dashed lines, respectively. 
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FIG. 3: (color online) The excluded production cross section 
for an SSM W' boson as a function of the W' boson mass. The 
excluded mass region is the region where the SSM W' cross 
section exceeds the excluded cross section. 
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FIG. 4: (color online) The excluded production cross section 
for an RS graviton as a function of the graviton mass. The 
excluded mass region is the region where the RS graviton cross 
section exceeds the excluded cross section. 



FIG. 2: (color online) The reconstructed WZ mass distri- 
bution for the data (points) with statistical error bars and 
the estimated Z+jets and diboson backgrounds, and the pre- 
dicted SSM W 1 distribution for M(W') = 600 GeV in the 
dilcpton channel. 



GeV. The W' signal acceptance and the expected back- 
ground yield are parametrized accordingly. A similar 
procedure is used to set limits on the RS graviton signal 
using the reconstructed transverse mass in single lepton 
events. We exclude an SSM W' boson in the mass re- 
gion between 180 and 690 GeV and an RS graviton in 
the mass range of 300 to 754 GeV at 95% CL. The gravi- 
ton mass exclusion assumes the dimensionless coupling 
parameter k/M p i =0.1, where k is the curvature scale 
of the warped extra dimension and M p i = M p i / v8~7r is 
the reduced Planck mass. The expected and observed ex- 



cluded cross sections as a function of the signal resonance 
mass for the combined single lepton, dilepton, and trilep- 
ton channels, along with the SSM pp — > W — > WZ pro- 
duction cross section are shown in Figure|3l The excluded 
cross section using only single lepton events and the RS 
graviton production cross section are shown in Fig. |4] 

In this analysis we assume a linear relationship between 
the resonance mass and total width, and that the width 
is smaller than the expected experimental mass resolu- 
tion. In some classes of models, the total width grows 
as a power of the mass yielding widths larger than the 
expected mass resolution. Using pythia W' —> WZ MC 
events with varying W' boson widths, we observe that 
our results are valid for widths below 10% of the reso- 
nance mass or, alternatively, for a coupling strength at 
the W'WZ vertex up to ten times the SSM value. 

In summary, with up to 5.4 fb -1 of Tevatron Run II 
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integrated luminosity we do not observe an excess of 
events over the SM background prediction for events 
in WW and W Z boson final states. We set limits on 
the production cross section multiplied by the branch- 
ing ratio for resonant WW and WZ boson pair pro- 
duction using two theoretical benchmark scenarios: SSM 
W' -> WZ and RS G -> WW production. Under these 
assumptions, we exclude an SSM W' boson with a mass 
between 180 and 690 GeV and an RS graviton with a 
mass between 300 and 754 GeV at 95% CL. Our novel 
use of the jet mass to select hadronic decays of the W and 
Z bosons was essential to obtaining such stringent limits, 
which are the best for these new physics scenarios. 
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